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ABSTRACT 

Context. Emission lines formed in decretion disks of Be stars often undergo long-term cyclic variations, especially in the violet-to-red 
(VI R) ratio of their primary components. The underlying structural and dynamical variations of the disks are only partly understood. 
From observations of the bright Be-shell star £ Tau, the possibly broadest and longest data set illustrating the prototype of this 
behaviour was compiled from our own and archival observations. It comprises optical and infrared spectra, broad-band polarimetry, 
and interferometric observations. 

Aims. The dense, long-time monitoring permits a better separation of repetitive and ephemeral variations. The broad wavelength 
coverage includes lines formed under different physical conditions, i.e. different locations in the disk, so that the dynamics can be 
probed throughout much of the disk. Polarimetry and interferometry constrain the spatial structure. All together, the objective is a 
better understand the dynamics and life cycle of decretion disks. 
Methods. Standard methods of data acquisition, reduction, and analysis were applied. 

Results. From 3 V/R cycles between 1997 and 2008, a mean cycle length in Ha of 1400 - 1430 days was derived. After each minimum 
in V/R, the shell absorption weakens and splits into two components, leading to 3 emission peaks. This phase may make the strongest 
contribution to the variability in cycle length. There is no obvious connection between the V/R cycle and the 133-day orbital period of 
the not otherwise detected companion. V/R curves of different lines are shifted in phase. Lines formed on average closer to the central 
star are ahead of the others. The shell absorption lines fall into 2 categories differing in line width, ionization/excitation potential, and 
variability of the equivalent width. They seem to form in separate regions of the disk, probably crossing the line of sight at different 
times. The interferometry has resolved the continuum and the line emission in Bry and HeI2.06. The phasing of the Bry emission 
shows that the photocenter of the line-emitting region lies within the plane of the disk but is offset from the continuum source. The 
plane of the disk is constant throughout the observed V/R cycles. The observations lay the foundation for the fully self-consistent, 
one-armed, disk-oscillation model developed in Paper II. 

Key words. Stars: circumstellar matter, emission line, Be - Stars: individual: ( Tau 



1. Introduction 

Circumstellar disks of classical Be-stars have been known for 
decades as the place of origin of their characteristic Balmer 
emission lines. Polarimetric ( |McLean & Brown] |1978| and 
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combined interferometric and spectro-polarimetric observations 
( IQuirrenbach et al.||T997] |Vakili et aL]|1998| l have brought the 
long discussion of the geometry of the disks to a definitive con- 
clusion. These observations confirm directly that the disks are 
geometrically thin with low opening angles of about 5-15 de- 
grees. The outer radius of the Ha emission was estimated at typ- 
ically 10-20 stellar radii. 

The detailed dynamics of Be-star disks, however, is not as 
clearly known. Various hypotheses exist, but the observations 
have not so far yielded fully unambiguous results. The relations 
between (a) stellar v sin i ([Hanuschik 1989) and (b) disk size 
( Quirrenbach et al.|1997[ ), and separation of the emission peaks, 
details of circumstellar line shapes (Hummel & Vrancken 2000 ), 
and the occurrence of central quasi-emission peaks in shell stars 
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(Rivinius et al. 1999) are all consistent with Kepler-like rotation, 
in which the rotation velocity varies as r~ J , where r is the dis- 
tance from the star. For a strictly Keplerian disk with circular 
orbits the exponent j is equal to 0.5. Analyses of spectral line 
profiles (Hummel & Vrancken 200 0} suggests that j < .65. 

|Hanuschik| ( |1994||2000| i and [Waters & Marlborough] ( fl994l l 
showed that any radially outward directed velocities in the disk 
must be lower than a few km/s. However, the angular momen- 
tum in a Keplerian disk increases with r 1 ^ 2 , Be disks are formed 
from matter outflowing from the star, and at least part of the disk 
eventually escapes from the star's gravity. Therefore, there must 
be a net angular momentum transfer from the star to the disk and 
then outwards through the disk. The mechanism for this momen- 
tum transfer, albeit not yet conclusively identified, is most likely 
related to dynamical viscosity. A viscous decretion disk model 
(|Lee et al.|1991| |Porter|1996} |Okazaki|200T]l has recently been 



successfully applied to the circumstellar disk of the Be star 6 Sco 
( Carciofi et al. |2006 1 and, at present, is the strongest candidate 
to describe the structure of Be-star disks. 

Purely circular Keplerian disks are also unable to explain the 
complex variability of the Balmer emission lines, notably the 
so-called V/R (the violet-to-red flux ratio of the emission peaks) 
variations. Most observations of V/R variations concern lower 
Balmer lines, but they can also be seen in virtually all other 
emission lines. Often, V/R variations are cyclic with timescales 
of 5 - 10 years (Okazak i|1991| l. The amplitudes can be large to 
spectacular. Nevertheless, many V/R cycles look strongly per- 
turbed because of concomitant large variations in the emission- 
line strength relative to the continuum and, presumably, the un- 
steady character of the mass loss from the central star. 

The first models for the V/R variations, a precessing ellipti- 
cal ring (Huan gjl973| ), spheroidal/ellipsoidal variable-mass loss 
(Doazan 1987 ), and a variable stellar wind causing expansions 
and contractions of the disk (Mennickent & Vogt 1991 ) were not 
physically self-consistent. They were finally made obsolete by 
Okazaki's model of one-armed density oscillations in the disk 
( |Okazaki|199T][T997l ). It is based on work by |Ka!^ ( [T983] l, who 
showed that the only possible global oscillation mode in an often 
thin Keplerian disk has m = 1. The idea of density waves was de- 
veloped further by|Papaloizou et aL| (| 1 992\ , [Okaz aki (2000), and 



Papal oizou & Savonije| ( |2006[ ). A few scattered interferometric 



observations detected decentered density enhancements in cir- 
cumstellar disks of £ Tau ( Vakili et al.|1998) > and y Cas (Berio 
|et al.|1999| ). 

The wide acceptance of the model is based mainly on qual- 
itative comparisons. Its detailed application to specific sets of 
V/R curves and emission line profiles has been pending to date. 
A strong reason is the poor repetitiveness of V/R cycles, which 
makes it difficult to decide what the model can be expected to 
reproduce. Furthermore, the typically long duration of the cy- 
cles has kept the statistics at a low level. Also, detailed compar- 
isons between model and observations have only become feasi- 
ble with the advent of three-dimensional radiative transfer codes 



( Carciofi & Bjorkman| |2006). The generic theoretical model- 
ing combining the dynamics of m=l disk oscillations with a 
radiative-transfer code was successfully pioneered by Hummel 
|& Hanuschikl ( p97l |. 

Stefl et al.| ( |2007[ ) analyzed the V/R variations in five Be-shell 



stars with a companion star. They find that in some stars the V/R 
variability is phase-locked to the orbital motions. In others it is 
not. A peculiarity seemingly limited to binary Be-shell stars (but 
not exhibited by all of them) is the appearance of triple-peak Ha 
emission profiles (which may also be described by a doubling 
of the self-absorption in the disk). £Tau is one of the examples 



presented by Stefl et al. (2007 1. The occurrence of such profiles 
seems restricted to a narrow phase interval at the beginning of 
the transition from V <k R to V ~ R. They are not reproduced by 
present versions of the disk-oscillation model. 

The star fTau (123 Tau, HR 1910, HD 37 202; B2IV) is 
among the most suitable targets for an in-depth observational 
and theoretical test of the disk-oscillation model. The present 
emission state, which started in 1990 ( |Guo et aLl [T995 their 
Fig. 2), shows very stable V/R variations with a cycle length of 
about 1500 days (Rivini us et al.|2006| ). Compared to other V/R- 
variable Be stars, secular and ephemeral variations have recently 
(see below) been small enough that observations obtained in dif- 
ferent cycles can be fairly safely assembled into one picture. 

For decades, ( Tau has been known as a spectroscopic binary. 
A comprehensive set of radial velocity measurements was com- 
piled and analyzed by Harmanec (1984 ), who derived a period 
of 132.9735 d. 

Being a bright object reachable from both northern and 
southern latitudes, £ Tau is one of the most observed Be 
stars. At the distance of 417 light years (corresponding to 
the Hipparcos parallax of 7.82 mas) not only the circumstellar 
disk can be resolved with present-day interferometers but also 
spectro-interferometry extends this resolution to the disk dy- 
namics. The limb-darkened photospheric diameter is estimated 
at 0.4 mas (see e.g. |Tycner et al.|2004]|Gies et al.|2007j ). At the 
Hipparcos distance, this corresponds to a radius of 5.5 to 6 R Q . 

The major axis of the Ho- emitting disk was measured by 
Quirrenbach et al. (| 1997 ) to 4.53 mas, corresponding to a radius 



of 11.3R+ and by Tycner et al. to 3. 14 mas (almost 8R+). Gies 



|et al.| ( |2007| ) derived the major axis of 1.99 mas (about 5 R*) for 
the K-band continuum emitting disk. These values were derived 
from the full width at half maximum (FWHM) of a Gaussian 
fit. The differences can be reconciled by assuming that the disk 
was less fully developed during the later observations as is also 
suggested by a lower Ha emission strength (see below). In the K- 
band continuum, |Gies et al.| measured the semi-major axis of the 
disk to 1.79 mas or 4.5 R*(Gaussian FWHM). All these values 
are lower than the estimated Roche-lobe radius for the primary 
of the system, which |Tycner et aL] compute as 5.3 mas (26.5 R*) 
at a binary separation of 9.2 mas (46 R*). 

This first paper of the series summarizes the extended obser- 
vational material (see Sect. [2]) and focuses on a detailed phe- 
nomenological description of the variability corresponding to 
the V/R cycle. The variations in the visible part of the spec- 
trum are analyzed in Sect. [3] The variability in emission lines 
in the JH TT-bands is described in Sect. |4] Sect. [5] is devoted to 
the VLTI/ AMBER spectro-interferometry and contemporaneous 
optical and IR spectroscopy, and Sect.|6]deals with new and pre- 
viously published polarimetric observations. The last two sec- 
tions discuss and summarize the observational results. 

The second paper (Carciofi et al., this volume) presents a de- 
tailed self-consistent physical model of the observations. It com- 
bines the 2D global oscillation model o f |Okazaki| (|1997|) and 
the 3D radiative-transfer code HDUST of ( |Carciofi & Bjorkman 
[20061 [20081 . 

2. Observations 

2.1. Visual Spectroscopy 

In order to maximize the coverage of the V/R Ha cycles, new 
observations were combined with data from the literature, the 
ESO Science Archive, and from the BeSS database (see below). 
Table [T] summarizes the spectroscopic data sets used for this 
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study. The descripti on of t he Heros and Feros spe ctrographs can 
be found in Kaufer ( 1998[ l and |Kaufer et al. ( T999}, respectively. 
The Flash instrument is an earlier version of the Heros spectro- 
graph with only one spectral arm (covering the range from 4050 
to 6700 A) and a somewhat lower resolution due to the use of a 
larger fiber, which determines the effective slit width. 

The OPD/LNA observation was made with a Cassegrain 
spectrograph, equipped with a 1200 grooves/mm grating blazed 
at 6562 A and a 1024x1024 pixel CCD. Observations at the 
Ritter observatory were done with an echelle spectrograph in 
the Cassegrain focus of the 1-m telescope, using the old cam- 
era equipped with a 1200x800 pixel CCD (cf. Sect.[2~4i. 

In recent years, the technological progress has brought 
bright stars like £Tau within reach of quite a few amateur 
telescopes and spectrographs. Spectra of many Be stars are 
regularly deposited in, and conveniently available from, BeSS 
(http://basebe.obspm.fr). As Table[T] shows, they account for a 
large fraction of the total optical spectroscopy used in this paper. 



2.2. Infrared Spectroscopy 

Fourteen infrared (IR) spectra were obtained between March 
5, 2004 and October 9, 2007 using the SpeX spectrograph at 
the 3.0-m NASA Infrared Telescope Facility (IRTF) on Mauna 
Kea ( |Rayner et £Ll[2006] l. With a C'3 x 15'.'0 slit in the short- 
wavelength cross-dispersed (SXD) observing mode, a resolving 
power of R~2000 spectroscopy was achieved from 0.8 to 2.4 
fim. The observing and data reduction techniques mirrored those 
used in previous SpeX programs to study classical Be stars (see, 
e.g., Wisniewski 2005 Wisniewsk i et al.|2007 1. Observations of 
£ Tau were immediately followed by observations of a nearby 
A0V star located at a similar air mass, to facilitate optimal tel- 
luric correction ( Vacca et al. 2003 ), along with a series of quartz- 
tungsten-halogen flat field and argon arc lamp exposures. The 
IDL-based Spextool software was used to perform the basic data 



reduction and spectral extraction using the techniques described 
~ 2004| l. 



in 



Cushing et al. 



2.3. Interferometry 



Our near-IR interferometric observations with AMBER ( |Petrov 
|et al. 12007) 1 and three 8-m telescopes (UT 1, 2, and 4) of ESO's 
VLTI interferometer were made during the night of December 
12, 2006. Interferometric fringes from all 3 pairs of telescopes 
(hereafter called baselines) were measured across the K-band, at 
medium spectral resolution of R w 1500. Because fringe track- 
ing was not supported at the time, it was not possible to properly 
record the full wavelength range. The observations were, there- 
fore, taken around the Bry (2. 16 pun) and Hei (2.06 /im) lines. 
The £ Tau observations were interlaced with observations of the 
calibration stars HD 39 699 and HD 59 686, for which diameters 
of 1.03 mas and 1.30 mas, respectively (Merand et al. 2005 [ l, 
were adopted. 

In principle, three distinct quantities can be extracted for 
each spectral channel of the AMBER spectro-interferometer: the 
fringe visibility (measuring the spatial extension of the emitting 
material), the fringe phase (measuring the position of emitting 
material), and the closure-phase (the sum of the 3 phases ob- 
tained for the 3 interferometric arms, sensitive to asymmetries 
in the source brightness distribution). In AMBER data sets, the 
signal-to-noise ratio on the closure-phase is very low and only 
visibilities and phases were used for the analysis. 

The data reduction was done using the amdlib-2.2 pack- 
age, which employs the P2VM algorithm (Tat ulli et al.| |2007 



In a first step, the instrumental calibration matrix with the stan- 
dard calibration files provided by ESO was computed. This is 
mandatory in order to properly convert the raw individual frames 
into raw interferometric visibilities and phases (about thousand 
measurements per observation). The next steps consist of aver- 
aging these different frames into individual measurements, and 
performing the final calibration. Based on our past experience 
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Fig. 1. Hff V/R variations (lower panel) and emission peak separation (upper panel)of £ Tau in the period 1992-2008 showing the 
detailed character of more than 3 V/R cycles since 1996. The symbols correspond to the following data sets: (+) Ondfejov slit 
spectrograph, (O) OHP 1.93m, (x) Pollmann, (H) Heros, (F) Feros, (a) Ritter Observatory, (■) BeSS database, (□) A. Carciofi, see 
also Table[T] Phases with triple-peaked emission line profiles are indicated by solid horizontal bars 



of spectro-interferometry with AMBER, two different strategies 
were adopted. In the following they are distinguished as absolute 
reduction and differential reduction. Only a brief description is 
included here, as both methods are now considered to be stan- 
dard data reduction steps for AMBER. 

Absolute data reduction: The objective of this method is to 
provide an absolutely calibrated estimation of the fringe visibil- 
ity (the fringe phase cannot be determined in an absolute man- 
ner). It requires the multiplicative effect of the atmospheric tur- 
bulence (generally called transfer function) to be properly esti- 
mated. 

- For each observation, the top 20% of the frames with the best 
signal-to-noise ratio were averaged. This ratio provides the 
most stable transfer function, see discussion in Tatull Fet al.| 
([2007} and |Millour etaLl ( |2007l ). 

- The transfer function was estimated by averaging all obser- 
vations of the calibration stars. Their scatter is a measure for 
the uncertainty of the transfer function. The resulting stabil- 
ity is very poor, with temporal fluctuations as large as 30%, 
partially due to the atmospheric turbulence and to the Unit 



Telescope infrastructures, which generate non-stationary vi- 
brations. 

- Absolutely calibrated quantities were derived by dividing - 
separately in all spectral channels - the raw visibilities of 
( Tau by the derived transfer function. Error bars on the 
transfer function completely dominate the final uncertainties. 

- Finally, all spectral channels were averaged (avoiding the 
He i and Bry lines) to provide a single visibility across the 
band for each observation. Formally it introduces a small 
amount of wavelength smearing, which is however negligi- 
ble at the level of precision. 



It is worth noting that this observation strategy has been opti- 
mized for differential interferometry, not for absolute calibration 
of the fringe visibility. Only a single calibration star has been ob- 
served for each AMBER instrument setup, and the observations 
of the calibration stars do not bracket properly the observations 
of the science target. Therefore it is hard to asses the reliabil- 
ity of the absolute calibration (usually done by comparing val- 
ues obtained using different calibration stars). Consequently, the 
fringe visibility derived with this method should be interpreted 
carefully. 
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Differential data reduction: This method aims mainly at a pre- 
cise measurement of the differential quantities (visibility and 
phase) in the Bry and Hei lines with respect to the adjacent 
continuum. Compared with the previous method, a significantly 
higher SNR can be obtained because the atmosphere affects all 
spectral channels in very nearly the same way. The main draw- 
back is that the continuum level is obviously lost. 

- First, as many as possible individual frames were averaged. 
To do so, all consecutive observations with the same instru- 
ment setup were merged together before averaging the 70% 
of the frames with the best signal-to-noise. It was tested that 
this method provides the best final accuracy on the differen- 
tial quantities. 

- The here presented AMBER data were corrupted by ripples 
with low beat frequency in the spectral direction (few cy- 
cles along the complete K-band spectrum) These artifacts 
mainly affect the phases and not the visibilities, explaining 
why they are not of a big concern regarding the absolute data 
reduction. To correct them, a high-pass filter along the spec- 
tral direction was applied to the phases. Different methods 
(high-pass filter, optimal filtering at the corrupted frequen- 
cies, manual fit of the continuum ripples) were tested but 
differences in the resulting He i and Bry phase profiles could 
be seen. 

- At that stage, the continuum levels for both the phase and 
the visibility are arbitrary. The brightness distribution of the 
K-band continuum emitting region has already been con- 
strained by |Gies et al.| l |2007) . Because they used longer base- 
lines and more precise absolute calibration, their results are 
significantly better than the AMBER measurements for the 
continuum (see previous paragraph). Therefore the AMBER 
continuum level around the Br y and He i lines was forced to 
to match the CHARA model of |Gies et al] ( |2007| l, which is 
displayed in Fig. 10 . 



2.4. Polarimetry 

The polarimetric data are from observations made with the 
Halfwave Polarimeter (HPOL; |Nordsieck & Harris |[T99cj] l. The 
instrument was used primarily on the 0.9-m telescope at the 
Pine Bluff Observatory (PBO), operated by the University of 
Wisconsin (UW). The data are available through the HPOL web 
site (http://www.sal.wisc.edu/HPOL) developed by M.R. Meade 
and B.L. Babler, and data prior to 1995 are also available through 
the NASA MAST archive. The database includes a wide range 
of hot stars, the polarimetric variability of which has been dis- 
cussed in Bjorkman & Meade (2005). 

Data from 1989-1994 were obtained using a dual Reticon 
array detector, which provided spectropolarimetry over a wave- 
length range of 3200-7600 A, with a spectral resolution of 25 A 
( |Wolff et al.|[l996l ). In 1995, the HPOL detector was upgraded 
to a 400 X 1200 pixel CCD camera, and data since then have 
covered a more extended wavelength range from 3200-10 500 A, 
with an improved spect ral resolution of 7 A below 6000 A and 



10 A above this point (Nordsieck & Harris 



1996). For further 



details about HP OL, see |Nook| (| 1 990| l and | Wolff et al] ([1996} 



As described by |Wolff et al. ( 1996[ ), the data were processed 
and analyzed using REDUCE, a specialized spectropolarimetric 



1 The interferences between the light beams reflected inside the 
AMBER polarizers created "Moire fringes" in the spectral direction. 
The optics responsible for these artifacts was identified and properly 
replaced in October 2008. 
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Fig. 2. V/R data from Fig. [T| phase binned with JD 50414 + 
1429 x E. Horizontal lines indicate the triple-peak phases (3-pk) 
from Table|2] which introduce a higher scatter in the correspond- 
ing phase interval because they make the definition of the V and 
R emission components doubtful. The letters in the upper part 
of the figure mark the phases of the AMBER observation and of 
repesentative visual and infrared spectra used in Sects. [33] and|4| 
(see also Tables [4] and [5} 



software package developed at UW. Instrumental polarization is 
monitored regularly as part of the observing program at PBO, 
and all data are fully calibrated for instrumental effects to an ac- 
curacy of 0.025% in the V band. 

For the purposes of this study of £ Tau, the spectropolarimet- 
ric data were binned to approximate broad-band UBVRI results. 
The individual uncertainty of an HPOL measurement was esti- 
mated by a statistics over a series of data of polarized standards: 
It is about 0.01 % for the polarization degree and 1 deg for the 
polarization angle in BVRI, and about two to three times as high 
in the f/-band. 



3. Visual spectral line variability 

3. 1 . Long-term evolution of the £ Tau circumstellar disk and 
its V/R variations 



|Harmanec| ( |1984| his Fig. 2) showed that the present well- 
pronounced V/R variations were absent in earlier decades: At 
least from 1920 to about 1950 there was no V/R cyclicity. A 
V//?-variable phase started in the mid-1950s and lasted for about 
three cycles until 1980. From 1980 to 1990 the star was again 
stable with V « R. From 1980 to 1985 the equivalent width 
(EW) of the Ha emission decreased from about -23 to -12 A 
while the star brightened by 0.3 to 0.4 magnitudes in the UBV 
passbands, also getting slightly bluer. This behavior is typical of 
Be-shell stars as described observationally by Harmanec] ( 1983) 



and predicted before by Poeckert & Marlborough 



(1978 



their 

Fig. 33) for the case of decreasing base density of the disk. 

From 1985 to 1990 the disk was in a low-density state ( |Guo| 
|et al.|[T995] their Figs. 2 and 3). The V/R variability was not 
well sampled observationally. But short-term (days) variability 
probably dominated (Guo et al. 1995 their Fig. 5). It would, 



then, have resembled a phase observed in n Cen by |Hanu schik 



et al. ( 1993), in which rapid V/R variations were due mainly to 



discrete mass loss events and the subsequent circularization of 
the ejected matter. 
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Fig. 3. Evolution of representative lines in the visible spectrum during the V/R cycle. See Table for the properties of spectra "A" 
to "F". The Ha V/R can bbe read off Fig. [2] 



Around 1990, the emission strength began to increase again Table 2. Properties of individual V/R cycles 

to a high level of almost -30 A, together with brightness and 

colour changes as expected for an edge-on disk with increas- 
ing base density. At the same time, the disk entered a phase 
of V/R variations. Since this replenishing of the disk, the emis- 
sion lines became weaker (Ha EW between -15 and -20 A ), 
but the V/R variations are still ongoing. The AAVSO database 
[^indicates that the magnitude was stable over the last 8 years 
(<m v > ~ 3. m 02 + 0. m 07). 



Cycle 


Max. V/R 


Cycle 


Triple peak 




number 


JD24 


V/R 


length [d] 


[JD 24 range] 


[d] 


I 


50414 


1.62 


1419 


51 137-51440 


303 


II 


51833 


1.58 


1527 


52580-53 094 


494 


III 


53 360 


1.44 


1230 


54005-54208 


203 


IV 


54 590 


1.40 


N/A 


N/A - N/A 


NA 



2 American Association of Variable Star Observers; 
htpp://www.aavso.org/data 
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3.2. Variations of the total Ha equivalent width 

Equivalent widths were measured by means of an automatic 
MIDAS ( |Grosb0l & Ponz|1990) procedure in the interval 6546 - 
6581 A. Except for a few outliers, all total EWs between 1992 
and 2008 are in the interval —11.5 to -20 A, with a mean of 
(-15.5 ±0.2) A. After an initial increase of emission lasting from 
about 1991 to 1994 the value remained rather stable. This is a 
further indicator of the validity of the assumption that the indi- 
vidual cycles are comparable to each other. 

In view of the uncertainties in the normalization of the 
echelle spectra with imperfectly corrected wiggles and system- 
atic differences between spectrographs with very different spec- 
tral resolution, a more detailed analysis of the measured to- 
tal equivalent widths and their temporal variations was not at- 
tempted. 



3.3. Shell absorption lines 

In the visual (and higher-resolution) spectra, two types of shell 
lines (dubbed "narrow-" and "broad-line group", NLG and BLG) 
and their associated cyclic behavior can be distinguished. Pure 
NLG-type shell absorption is found in lines of Fen (near-UV 
lines only), Fein, Nin, Oi, Nai, and Hei. The cores of the 
Balmer lines, too, fall into the NLG category. Pure BLG-type 
shell absorption is seen in some lines of Si n, Fe n (visual wave- 
lengths), Cm, Tin, Mgi, and Oi. This group also includes 
Mg ii 4481, which is among the easiest lines to excite in a Be 
circumstellar environment. The NLG/BLG classification of the 
most important spectral lines in the optical spectrum is marked 
in Tabled] 

The assignment of a line to BLG or NLG is not in all 
cases and not at all epochs unique. For instance, Fen5137 and 
Si ii 6347 normally show the BLG characteristics, but may have 
additional NLG components in certain phases of the V/R cycle. 
Other BLG lines, too, exhibit additional temporary NLG com- 
ponents, especially around V/R phase t w 0.25 (see Sect. 3.5 1. 



Table 3. Selected spectral lines in the optical spectrum and 
their classification into broad-line (BLG) and narrow-line group 
(NLG), respectively, as described in Sect. 



3.5 



Typical representatives are the strongest shell lines due to Si n, 
Ca n, Ni n, and Fe n. 

Besides the obvious property of their widths, the two groups 
also differ in their variations of radial velocity and symmetry as 
described in Sect. 17. 71 



Wavelength Ion Comment 



a i on ha /z 
j 1 o / . /40 


T4o t 

rie I 


iNLLr, no pnot. 


Q/l/^C A/IO 

j40j.04Z 


Mi ii 
INI II 


MT C 
INLLr 


Q/1AQ £7Q 

j40o.O/ o 


Fe ii 


MT C 
INLLr 


J4 1 1 JOO 


Mi ti 
INI II 


MT C 
INLLr 


ocn OC7 

jjl j.yo / 


JNl II 


NLG, strong 




rie I 


pure phot. 


1 A 1 O 

jjj4.41Z 
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pure phot. 


•5C07 OAG 

jjo /.Zoo 


T4o T 

rie I 


pure phot. 


JOj4.ZjO 


rie i 


pure phot. 




r^o ti 

La II 


DT (• . "VTT 

dIAj + INLLj 




Psl II 


pure phot. 


/inno ^oo 
4UUZ. jyZ 


ol II 


DT (' 

dLu 


A AA7 (YX 1 

4Uo / .UjI 


JNl II 


di r i mt 
13 Lu + INLLj 


41Zo.Uj4 


Ol II 


RT (' _i_ MT C* 


A 1 1C\ 071 

41 jU.o /Z 


c; ti 

Ol II 


DI P i MT f~* 

rSLLr + INLLj 


/lO^Q 1 ^/1 
4ZJo. 1 J4 


Fe ii 


DT {• 

dLu 


4zy4.uyy 


t; tt 
ll II 


D T f • 
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A A C 1 1 OA 
44 o 1 . 1ZO 


Mg ii 
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dLu 


ARIA 1 97 

H-oZH-. 1Z / 


V I 11 


RT C 
DLU 


cm K A7£ 
JU1 J.O / o 


rie i 


inllt + weaK pnoi. 


c 1 nn TV7 
J 1UU. /Z / 


Fe ii 


DT f ' 

13 Lu 


J 1 Z /.jo / 


Fe in 


MT C 
INLLr 


C 1 CA 1 1 1 
J 1 JO. Ill 


Fe in 


MT f ' 

INLLr 


c 1 AO CYX'X 
jlDV.UjJ 


Fe ii 


RT (' _i_ MT C 
13LLr + INLLj 


c 1 70 AC/1 
J 1 /Z.0o4 


Mgi 


RT C 
13 Lu 


cioq AHA 
J 1 OJ.DU4 


Mgi 


RT C. 
13 Lu 


c 1 07 C77 

jiy /.j / / 


Fe ii 


RT (' _i_ MT C* 
13 Lu + INLu 


^77A ^Al 
JZZO. J4J 


Ti TT 

ii ii 


RT C 
D 1 A J 


K'X 1 A A1 ^ 
J J 10.01 J 


Fe ii 


RT C. _i_ MT CI 
D Lu + INLO 


5679 558 


Nn 


pure phot. 


5875.625 


He i 


NLG + phot. 


5889.951 


Nai 


NLG 


5895.924 


Nai 


NLG 


6147.741 


Fe ii 


BLG 


6149.258 


Fe ii 


BLG 


6156.755 


Oi 


BLG 


6158.187 


Oi 


BLG 


6347.109 


Sin 


BLG + NLG 


6371.371 


Sin 


BLG + NLG 


6678.154 


He i 


NLG + phot. 


7774.166 


Oi 


NLG 



3.4. Length of the V/R cycles 

Some confusion may occur due to different definitions of the 
V/R ratio in the literature. They differ in the reference flux level, 
with respect to which the peak maxima are measured. This pa- 
per uses the more common definition (V/R = Fv/Fr, where 
Fy and Fr are the relative fluxes at the violet and red emis- 
sion peak maxima, respectively), according to which V/R is 
computed without prior correction for the underlying continuum 
level. Therefore, the results are less influenced by the spectrum 
normalization than in the other definition, which uses only the 
flux above the continuum level. Consequently, this approach en- 
ables to combine reasonably well V/R values even in hetero- 
geneous databases. Since the photometric variability has been 
small (Sect. |3.1 1, while the H a line emission was strong, the 
continuum variability should have little effect on the V/R values 
measured in this way. The V/R values obtained for the interval 
1993-2008 are visualized in Fig.[T] 

The observations before JD 2 450 500 are insufficient to con- 
clude whether the V/R cycle was not yet stabilized or whether 
the sampling was too sparse to deduce a more regular V/R vari- 



ability. However, the subsequent observations define nicely three 
complete cycles marked I, II, and III in Fig.[T] Cycle IV just com- 
menced about JD 2 454 590. Starting Julian Dates and lengths of 
the cycles are summarized in Table [2] The cycle-averaged V/R 
curve is displayed in Fig. [2] 

The mean V/R curve is smooth and roughly symmetric ex- 
cept for some time shortly after the minimum. During this as- 
cending branch, the Ha emission is split into three peaks (or the 
self-absorption is split into 2 components), making the definition 
of 'the' V and 'the' R component ambiguous (see Sect. 3.5.4 
and Figj6]». An increased scatter in V/R values is the result (see 
Figs.[T]and|2]). Obviously, the disk variability cannot be charac- 
terized by a single parameter such as V/R alone. 

By contrast, the V/R maxima are very well defined and can 
be used as fiducial marks for measuring the length of V/R cy- 
cles. Fourth-order polynomials fitted to the peaks over different 
baselines indicate an error in their position of 5-8 days. 
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In spite of the systematic uncertainty in V/R, Fig. [T] seems 
to suggests a secondary V/R maximum that coincides with the 
triple-peak phases and follows the main minimum. Apparently, 
such secondary maxima in the cyclic V/R variations were not so 
far reported for any Be star. Their height and widths vary from 
cycle to cycle and the strongest secondary peak can be recog- 
nized in Cycle II at JD 2 452 700 - 2 453 000. 

A formal Fourier frequency analysis was only applied to 
Ha V/R values not affected by a third component (see Table|2] 
for the corresponding dates). Such a restriction is suggested by 
the sensitivity of most period-search methods to non-sinusoidal 
perturbations and to a lower accuracy of the V/R values dur- 
ing the triple-peak phase. The CLEAN method of the sinusoid 
fitting ( |Kaufer et al]|1996| l gives periods P x = (685 + 5) d 
and P2 = (1405 + 26) d, obviously the former one is the har- 
monic produced due to variable cycle lengths. Scargle's method 
(Scargle 1982) is less likely to return harmonics and suggests 
1429 days. 

This mean duration of the V/R cycles is shorter than the val- 
ues derived by Pollmann & Rivinius (2008 1, 1475 days, and by 
Rivinius et al. ( |2006[ l, 1503 days. The differences are well ex- 
plained by the different data sets used. In particular, the earlier 
studies could not include the shorter-than-average triple-peak 
phase of Cycle III. Even more different and variable cycle length 
can be identified during the 1955 - 1980 V/R variable phase. The 



Spectrum E 



compiled data presented by |Harmanec ( 1984] > show three max- 
ima, defining -in the same sense as in this paper- two complete 
cycles of about 2290 and 1290 days. 

Radial-velocity (RV) measurements suffer from the inhomo- 
geneity of the database in a similar fashion as the absolute emis- 
sion strength does. The impact of this inhomogeneity is reduced 
in the RV separation of the red and violet emission peaks (Fig.[T| 
The separation follows the V/R cycle, but with a possible phase 
shift and a considerably higher scatter, which not surprisingly is 
most prominent when the triple-peak profiles are present. The 
dates of maxima of peak separation are delayed by 100 - 200 
days with respect to the V/R maxima. 
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Fig. 4. Ha, H/3 (solid), and Oi8446 (dashed) line profiles of 
spectrum "E" (left) and the AMBER-simultaneous FEROS spec- 
trum during triple-peak phases. The absorption on the blue side 
of 1 8446 is a Paschen line 



ably optically thin lines like Fe 11 63 18 has just finished in spec- 
trum "F" and is about to return after spectrum "E". (Note that 
in Fig. [2] "F" appears at an earlier phase than "E". But this is 
only an artifact of the co-phasing of cycles with different length. 
Phenomenologically, "F" is actually following "E".) 

This classification is restricted to the visual spectral range. 
For the IR spectra, which are more limited in phase coverage 
and spectral resolution, a separate scheme is introduced in Sect. 

a 

The following subsections describe the main optical V/R 
phases and their transitions. The definition of the V/R phases is 
for Ha. The variability of other lines may have significant offsets 
in phase, as detailed below. 



3.5. Spectral evolution during the V/R cycle 

The repeatability of the V/R cycles for more than a decade justi- 
fies the introduction of the concept of an ephemeris and a phase 
so that observations can be compared between cycles. In the fol- 
lowing, all phases refer to the formal cycle length of 1429 days 
from above and adopt the V/R maximum of JD 2450414 as the 
reference epoch. The co-phased V/R curve is given in Fig. [2] 

In order to more completely characterize the disk variabil- 
ity than the V/R ratio can do, 6 optical reference spectra were 
selected. Their V/R phases are marked in Fig. [2] The spectra 
were sorted such that the phase of equal peak height in presum- 



Table 4. Optical spectra representative of the V/R-cycle. Flash 
is the predecessor instrument of HEROS, i.e. the same instru- 
ment before the extension with a blue channel 



Spectrum 


[JD24...] 


Instrument 


Comment 


A 


53 333.8 


UVES/ESO Archive 


nightly average 


B 


52009.3 


HEROs/Data Set C 


single spectrum 


C 


53 700.8 


UVES/ESO Archive 


night average 


D 


52 364.3 


HEROs/Data Set C 


single spectrum 


E 


52 725.3 


HEROs/Data Set C 


single spectrum 


F 


48 347.6 


Flash/LSW Archive 


single spectrum 



3.5.1. V >R 

This phase is represented by the "A" profiles in Fig. [5] In all 
emission lines the violet peak is higher than the red one. The 
shell absorption lines weaken and become narrower. Both NLG 
and BLG are red-shifted with respect to the systemic velocity. 
However, there are differences: The red edge of BLG shell ab- 
sorption is very steep, and the blue edge rather asymptotically 
joins the continuum. In the NLG, it is the blue edge which is 
steeper than the red one. All ionic species are present but some 
low-ionization lines begin to weaken. 

3.5.2. V = R, descending to V < R, and deep absorption 

When the strength of the R peak in Ha begins to approach the 
one of the V peak, the peaks in H/3 are already of almost equal 
height, and Fe n lines even exhibit an inverted V/R ratio (cf. spec- 
tra "B" and "C" in Fig . [3j . 

Some of the low-ionization shell lines of the BLG category 
have completely vanished, while the NLG shell lines are at max- 
imum depth now. The Balmer series is maximally visible up to 
H40/41. At other phases, the limit is reached around H34/H35 
(see Fig. [7]). 

While pure BLG lines such as Mg 11 448 1 are slightly blue- 
shifted, pure NLG lines reside at zero velocity and are very nar- 
row. Also shortwards of the Balmer discontinuity several shell 
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Fig. 6. Evolution of the Ha triple-peak profiles in the two most 
recent V/R cycles of ( Tau. The shifts between the spectra cor- 
respond to differences in phase as indicated on the left side. 
The phase is computed separately for each cycle, using differ- 
ent starting dates (V/R maxima) and lengths as listed in Table [2] 
The vertical bar indicates the 100 % continuum level 



lines are far more clearly present than at other phases, phe- 
nomenologically belonging to the NLG. Their FWHM is only 
about 17 to 20km s _1 , i.e. about the thermal width. 

Shell lines of Fen and Sin, however, develop a two- 
component structure related to both groups: a very narrow and 
deep absorption core at zero velocity is superimposed on a broad 
and shallower absorption that has its own distinct core farther to 
the blue. The two cores can be clearly identified as parts of the 
BLG and NLG characteristics, respectively. 

3.5.3. V <R 

The "D" profiles in Fig. [3] illustrate this phase. The red peak 
reaches its maximal strength. Both NLG and BLG lines are 
blue-displaced with respect to the systemic velocity. The dual 
NLG/BLG characteristics of some lines weakens and finally 
vanishes in that NLG cores disappear while the remainder of 
the profile evolves as if belonging to the BLG. 

3.5.4. V = R, ascending to V > R, and triple-peak structure 

Spectra "E" and "F" in Fig. [3] offer a general overview of this 
phase, and Fig. [6] illustrates the temporal evolution in Ha. The 
most intriguing property of this phase is that the Ha line emis- 
sion strongly deviates from the classical Be star profile in that 



it no longer shows clear double peaks with a well-pronounced 
central depression. Rather, the central depression is filled in (or 
split) and three peaks of similar strength co-exist. Although the 
central depression does not get filled in in the other Balmer lines, 
and they do not show multiple peaks, their structure is more com- 
plex than in other phases. 

The 1 8446 line, which is excited by fluorescence of Ly/3 
and therefore usually assumed to present an optically thin tracer 
of the hydrogen distribution, does not show any sign of a triple 
peak (see Fig. HJ. Instead, the profile is double-peaked with a 
central depression. Noteworthy is that the radial velocity of the 
red peak of the O i profile coincides with the local minimum in 
the triple peaked Ha and H/3 profiles. 

The shell absorption lines are strong and yet relatively broad. 
Lowly ionized species like Ti n, Cr a, or Mg n prevail, but high- 
ionization species like Fe in still produce noticeable absorption. 

Spectrum "F" (Fig. [3j was obtained in 1991. Since then, 
the shell spectrum has generally weakened. But there is no 
doubt that spectrum "F" is representative of the general V/R pat- 
tern about this phase. The latter includes that, between spectra 
"E" and "F", BLG shell absorptions shift from zero to moder- 
ately positive velocities. Examples are Mg n4481, Si n6347, and 
Fen 53 17 in Fig. [3] This is accompanied by a change in profile 
symmetry. While in "E" there is a sharp blue edge, "F" exhibits 
a sharp red edge. 
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Table 5. Representative infrared spectra 
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Wavelength [A] 



1.15 



1.05 



0.95 




3465 



3470 

Wavelength [A] 



3475 



Fig. 7. The region of the Balmer discontinuity (upper panel) and 
the narrow shell lines (lower panel) in the Balmer continuum in 
reference spectra "A" (lower plots) and "C" (upper plots), see 
Sect. Ed 



In the NLG shell lines, however, the position of the line min- 
imum as well as the symmetry change in exactly the opposite 
way, from low positive velocity to zero velocity in the same two 
observations. See, e.g., Fera5156 and He 1 6678 in Fig. [3] 

There are very few pure emission lines, e.g. Fen6319. But 
all of them, as well as lines such as H/3, where emission still 
dominates over the shell absorption, already show the V peak 
higher than R. That is, the V/R phase of these lines is ahead of 
the one of Ha. 

Fig.[6]illustrates the evolution of the triple peak profiles in the 
two most recent V/R cycles. The triple-peak interval in Cycle II 
lasted for almost 500 days, while in Cycle III it was only 200 
days. As can be seen from Figs. [T] and [2] the onset of the triple- 
peak interval is typically fairly close to phase 0.5. Much stronger 
cycle-to-cycle variations are evident for its terminal phase. In 
spite of the varying time scale, however, the morphological evo- 
lution of the Ha line profile follows approximately the same pat- 
tern in all cycles. 



4. Infrared spectral variations 

As for the optical domain, a number of phase-representative ref- 
erence IR spectra were selected to facilitate the description of 



Spectrum 


Date [JD 24...] 


Comment 


I 


53 722 


22 days after spectrum "D" 


II 


53 982 


just before triple peak 


III 


54097 


contemporaneous to AMBER 


IV 


54281 


just after triple peak 


V 


54 382 





the main variability. Their dates are marked as I through V in 
Fig. [2] and listed in Table[5] 

The IR spectra (see Fig. [8} cover the JHK bands and the Ha 
V/R phases from about 0.3 through 0.8. Although of lower res- 
olution, some of the conclusions already drawn from the visual 
spectra are confirmed: The V — R deep absorption phase has 
the same characteristics in the IR, as far as H i and He i are con- 
cerned. In the triple-peak phase, the shell absorption gets shal- 
lower in He i and even disappears completely in the higher-order 
Brackett lines, leaving a pure emission profile although with a 
pronounced central depression. All observed Brackett lines are 
double-peaked at all times. But Bry and Br£ develop single-peak 
profiles when the one in Ha shows a triple-peak structure. Also 
metallic lines, in particular those of Fe n may show single peak 
profiles, but at other phases. But it must be suspected that this 
merger is only apparent due to the limited spectral resolving 
power. 

The IR lines provide the basis for a more in-depth discussion 
of the phase lags between Ha and other lines (Sect. 3.5.4 1. At 



some phases, they are well visible even within Hi emission of 
the Brackett series. In Spectrum III of Fig.[8]the lower Brackett 
lines still show V < R and, in Bry, even V <K R, but towards the 
end of the series V/R reaches unity. In a spectrum taken less 
than a month later, some of the latter lines even appear with 
V > R while Bry looks unchanged. The V — R state is also 
reached in He i, Fe a, and Mg n. These two particular observa- 
tions (incl. Spectrum III) were obtained almost simultaneously 
with the AMBER data described below. 

The phase lag between higher- and lower-order Brackett 
lines is also well visible in Fig. [5] Already the first IR obser- 
vation on JD 2 453 069 shows that, while log(V/R) is still below 
unity for Ha, it is about unity for Bry and Br5, and still larger 
for higher Brackett lines. In general, the V/R values for both Bry 
and Br<5 precede the V/R curve of Ha, but clearly lag behind the 
ones of the Br 12-16 lines. The V/R amplitude of the Br 12-16 
lines is slightly lower than the ones of Bry and BiS. But it is 
lower by a factor 8-10 relative to the one of Ha. 

Contrary to Ha (see Fig. E|, a distortion of the V/R curve 
around phase 0.6 due to triple-peak profiles cannot be recognized 
in any IR lines. Observations at higher spectral resolution are 
desirable to establish this definitively. 

A very conspicuous type of behavior without counterpart in 
the visible range is associated with the infrared Ci lines (see 
upper right panel of Fig . [8]> . IR multiplets 1 and 24 of this species 
are seen in emission in six out of fourteen IR spectra (JD24: 
53 069, 53 639, 53 722, 53 723, 54038, and 54 070) and absent 
in the others. Since intermediate cases were not observed, this 
may constitute some type of "on/off" behavior. Contrary to the 
typical circumstellar lines described above, these lines have a 
symmetric single-peak profile and exhibit no variability other 
than their being either present or absent. No relation could be 
found of the presence or absence of the C i emission to the V/R 
cycle or the 132.9735-d orbital period. It was also checked that 
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Fig. 9. Visual spectroscopic state of ( Tau during the AMBER observations 



the variations cannot be produced by an overlapping OH telluric 5. Interferometric results 

emission (Rous selot et al.|2000) . _ u , .^.n^r, 

1 1 ' 5. 1 . Results from AMBER 



In Fig. 10 the left panel shows the coverage of the uv plane by 
our AMBER observations. Following the ESO/VLTI documen- 
tation, the following baseline position angles were accepted: Ul- 
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U3 = 32.369, Ul-U4 = 60.396 and U3-U4= 110.803 deg. The 
right panel shows the absolutely calibrated squared visibilities 
in continuum. Because the transfer function is very uncertain, 
the final precision is rather poor, errors are as high as 20-30%. 
The model visibilities for baselines along the major and minor 
axis of the ( Tau model derived from the CHARA K'-band ob- 
servations by |Gies et aL| ( |2007| l are over- plotted. The AMBER 
continuum observations - in spite of the large scatter - are com- 
patible with the CHARA model. However their accuracy and the 
fact that £ Tau is only marginally resolved by AMBER/VLTI in 
the K-band continuum prevent us to test the model in more de- 
tail. 

AMBER spectro-interferometric measurements were com- 
bined with the CHARA model in order to derive visibilities and 
phases in Bry and He i 2.06 // emission lines (see Sec. 2.3 1. The 
differential values are shown in Figs. 1 1 and Figs. 12 for the re- 
gions around Bry and Hei 2.06//, respectively. In all baselines, 
the signal is much stronger in the Bry line (2.17 //m), but it is also 
clearly detected in the He i line (2.06 //m). The drops in visibility 
prove that the Bry and He i emitting regions are more extended 
than the continuum emitting region. Moreover, an asymmetry 
in the visibility profile is clearly seen, with the visibility in the 
red part being lower than in the blue part. Fitting these visibility 
profiles across the lines requires proper modeling. Such a work 
is presented in paper II. 

Apart from the visibility signature, also a phase effect was 
found in all baselines for Bry, indicating that the photocen- 
ter across the line is shifted with respect to the photocenter 
in the continuum. To allow a geometrical representation of the 
AMBER phases, the differential -phases <p(X) across the Bry line 
were converted into 2D astrometric shifts p(A) by inverting the 
well-known formula for marginally resolved interferometric ob- 
servations (Lachaume 2003): 

m = -2n-p{X)-~ (1) 

A 

where b is the interferometric baseline vector projected onto the 
sky. Uncertainties were propagated to the astrometric vector p by 
standard formulas. Neither the uncertainty in the baseline length 
nor in the spectral calibration were taken into account, because 
they affect all spectral bins in the same way. 

Figs. 1 1 and 12 prove the presence of phase differences of 



400 - 



200 - 



the V and R emission components with respect to the contin- 
uum. The observed (differential) phases are the flux weighted 
vector sum of the source phases of the continuum and line com- 
ponents. With a physical model of the disk and radiative trans- 
fer, the fluxes are known and the phase differences between the 
continuum and line emitting regions can be translated into their 
angular separation. Such a model is presented in Paper II. The 
model-independent conclusion is that the line emitting regions 
and the continuum source lie in one common plane but that the 
V and the R components arise from different locations that do 
not coincide with the continuum source. 
Fig. 
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left panel shows the computed relative offsets in 
angular units. The phases across the Bry line is converted into 
2D photocenter shift in the plane of the sky. The photocenter is 
displaced toward the NW direction in the blue part of the line, 
and toward the SE direction in the red part of the line. This is 
the clear signature of the rotating disk emitting the Bry line. The 
position angle of the displacement is perfectly compatible with 
the position angle of the CHARA model based on measurements 
made in the K' continuum, as shown in the right panel of Fig. 1 3 
Moreover, the amplitude of the photocenter displacement is sig- 
nificantly larger in the red wing (280 //as, SE) than in the blue 
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Fig. 10. Schematic comparison of the VLTI/AMBER observa- 
tions with the CHARA model by |Gies et aL"1 ( [2007| l. In the colour 
plot available in the on-line version of the paper, the green, red 
and blue colours correspond to the UT3-4. UT 1-3 and UT 1-4 
baselines, respectively. Left panel: Coverage of the uv plane by 
AMBER observations of Dec 12, 2006. The sub-frame in the up- 
per right corner shows the CHARA model projected on the sky 
(before the transformation). The sub-frame covers the region of 
about 3mas x 3mas and corresponds to the K-band continuum. 
The squired visibility of the AMBER observations is linearly 
proportional to the size of symbols, for the CHARA model (in 
the sub-frame and transformed in the background) to the bright- 
ness. Right panel: Absolutely calibrated continuum squared vis- 
ibilities as a function of the baseline length. Squared visibili- 
ties of the major and minor disk axes derived from the CHARA 
model are overplotted. 



wing (120 //as, NW). The simplest explanation is an asymme- 
try in the Bry-emitting material, with the SE part being brighter 
and/or more extended than the NW part. Again, modeling of this 
signal in the framework of the one-armed disk-oscillation model 
is developed in Paper II. 



5.2. Optical and IR spectra accompanying the AMBER 
observations 

When AMBER observed £ Tau, simultaneous FEROS optical 
spectra as well as quasi-simultaneous IRTF spectroscopy were 
also obtained. Fig. [9] depicts representative visual line profiles. 
Further FEROS observations were made throughout December 
2006, which differ but little from the one at the time of the VLTI 
observations. In the IR, Spectrum III (see Fig. [8]) comes closest 
(12 days before) in time to the AMBER observations. A next 
spectrum, obtained 16 days after AMBER data, is very similar. 
Therefore, the low variability in both optical and IR spectra sug- 
gests that the selected spectra represent well the dynamics of the 
disk at the time of the AMBER observations. 

As IR spectrum III shows (Fig. [8]), the emission is almost 
symmetrical in He 1 2.06 //m, so that the much lower amplitude 
of the phase signature in the same line is in good agreement with 
the spectroscopic data. The distribution of the He 1 2.06 //m emit- 
ting material is symmetric with respect to both systemic velocity 
and location of the continuum source. 

The situation is different for Bry but the reasoning is the 
same: The closure phase is undoubtedly non-zero for this line 
and proves a spatial asymmetry in the brightness distribution of 
the system (continuum vs. line-emitting region). At the same 
time, the V/R value of Bry is far from unity. The combination 
of these two asymmetries places the bulk of the Bry-emitting 
region to the receding sector of the disk. 
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Fig. 13. Left: Photocenter shifts derived from the AMBER rel- 
ative phases across Bry. The maximum shift is about 0.4 mas 
within the plane of the circumstellar disk, while no significant 
offset perpendicular to it can be found (black line). Right: The 
position angle derived from our differential data overplotted on 
the model of |Gies et al.] ( |20"07] >. 
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The scheduling of the AMBER observations was very suc- 
cessful in that they could be made during a triple-peak phase of 
the Ha emission. However, with single-epoch data it is difficult 
to say whether they harbor any related hidden anomaly. For a 
discussion of the azimuth (in a not point-symmetric model) of 



the region of formation of this spectral structure see Sect. 7.6 



6. Polarimetric results 

Since the V/R cyclicity before JD 2 490 000 is not well con- 
strained by the observations and, as described in Sect. 3.2 £Tau 



reached a stable disk state only after that date, only polarization 
measurements obtained after JD 2 490 000 were taken into ac- 
count for the computation of the mean values shown in Table [6] 
The polarization angle (PA) is the same in all bands within 
both the individual uncertainty and the standard deviation. 
Although the standard deviation of the PA measurements is 
slightly higher than the uncertainty estimated for a single mea- 
surement, this difference is not large enough to conclude a vari- 
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able polarization angle. Even after inclusion of all data, i.e. also 
the observations taken before the V/R cycle is well defined, the 
values for the PA in Table [6] would change only in the last digit, 
if at all, so that the conclusion about a non-changing angle is 
valid over the entire period of observation. 

Other than the angle, however, the polarization degree does 
vary. The variability shown in Fig. [14] can be disentangled into 
two components. A slow and secular change is visible well in 
the VRI and, though somewhat noisy, also in the U B bands in 
Fig. 



14 This behaviour is in agreement with the general evo- 
lution of the Ha disk as described in Sect. |3.2| The mean po- 
larization degree is lower before JD 2 448 000, then rises until 
JD 2 450 000, and after that remains stable again. Superimposed 
on that pattern is short-term variability. Unfortunately, the data 
do not cover the timescales required to properly constrain the 
short-term component. However, the observed behavior is not 
in contradiction with polarimetric variability produced by a se- 
ries of discrete and individually rather minor ejections of matter 
into the immediate stellar environment, as already suggested in 
Sect. ED 

In summary, neither the variability of the polarization degree 
nor the polarization angle, which is stable, can be connected to 
the cyclic V/R variability. No modulation due to the orbital pe- 
riod could be detected. 



7. Synopsis and discussion 

7.1. ( Tau spectral type 

Table [3] lists several near-UV He i lines as having purely pho- 
tospheric profiles. They offer the opportunity to shed some light 
on the various discrepant spectral types that have been published 
for £Tau. A comparison of the high-quality UVES data to sev- 
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eral early B-type stars in the UVES POP-library (Bagnulo et al. 
|2003] ) showed unambiguously that the observed strength of these 
lines is inconsistent with an effective temperature corresponding 
to a spectral type of £Tau later than B2. However, an exact re- 
determination of the spectral type is beyond the scope of this 
study. 

7.2. Binarity 

With f Tau having an unseen companion, any variability must 
be checked for a relation to the 132.9735 d orbital period estab- 
lished by Harmanec ( 1984| l. However, the strength of the disk 
during the observations implies that the vast majority of the spec- 
tral lines are contaminated by line emission primarily varying on 
a different time scale. The other lines are mostly too shallow for 
radial-velocity measurements that would permit a verification of 
Harmanec s amplitude of lOkm/s to be attempted, which 
was derived from observations at more V//?-quiescent epochs. 
The photospheric He i lines are all in the blue spectral region, 
which is covered only by a small number of spectra used in this 
study. 

It follows that the used dataset of optical spectra, focused 
on the Ha V/R variations is not suitable for a refinement of or- 
bital parameters. As a compromise only a check of the period 
using the Hei4026 line was performed. Radial velocities were 
measured by fitting Gaussians to the line core. Phasing this data 
with [Harmanec 's ephemeris gives a plot similar to his Fig. la 
but with higher errors and a few outliers. A time-series anal- 
ysis finds the second strongest peak corresponding to a period 
of 131.6 ± 0.5 d, in reasonable agreement with Har manec] The 



strongest peak present is a nearby alias of this one, probably en- 
hanced by noise effects. No other quantity investigated in this 
work shows a significant modulation with the orbital period. 

Even over the given large spectral range, lines due to a sec- 
ondary, either hot or cool, were not found. In particular, no 
trace of He n 4686 was detected that would have indicated a hot 
subdwarf-like companion. 



Re-analysis of the publicly available photometric data ( Guo 
let al.|1995[|Pavlovski et al.|19971|Harmanec et al.|1997| and data 
from the HIPPARCOS-mission) did not confirm the presence of 
photometric eclipses suggested by |Bozic & Pavlovski| p988 ). 
Rather, the fainting at phase zero they observed seems to be an 
event unrelated to the orbital phase. There are more such fadings 
but at random phases and none of them at phase zero. 



7.3. Disk orientation on the sky 

In an axi-symmetric circumstellar disk the PA is in general per- 
pendicular to the disk plan^jand thus parallel to the disk position 
angle (x). 

The spectropolarimetric monitoring of ( Tau from the PBO 
observatory shows that the polarization angle has been remark- 
ably constant from at least 1989 to 2004 (see Sect. 15). The av- 
erage V-band polarization angle is PA PB o = 31.3 + 1.5 ° with 
an individual measurement uncertainty of 1 deg. The data from 
Limber Observatory, taken between 1984 and 1997 (McDavid 
1 1999 1, gave very similar results: PA L i m b er = 32.5 + 1.1°, and a 
typical individual uncertainty of 0.5°. 

Further, there are three independent measurements of the 



position angle from interferometric studies. Quirrenbach et al. 



( |1997| l determined from a 2D Gaussian fit for data from the Mark 
III interferometer a value of ;f Markin = 31 ± 4° . Using data from 
NPOI |Tycner et al.| ( [20"04"] > obtained *npoi = 28 ±4°, while [Gies 
[etiLl ( [2007) report |y C HARA = 37 + 2° from their CHARA data. 
Finally, the VLTI data of Sect.|5]show that Camber = 32 ± 5° 

The orientation of the disk on the sky, as indicated by all 
these measurements, has not changed since 1984, and in particu- 
lar there seems to be no binary phase-locked variability through- 
out the V/R cycle, at least within the individual measurement 
uncertainties. The stability of the disk position angle is typical 
for Be stars, the only few counter-examples ever found are sus- 



pected to undergo precession of a non-equatorial disk (Hummel 
|T^|Hlratal2007l ). 

7.4. V/R cycle length 

The analysis of the up to now most comprehensive observational 
data set of ( Tau confirmed that its V/R variations follow a cycle 
with relatively stable amplitude and length of 1405 - 1430 days 
duration during the present V/R variable phase starting at the 
beginning of nineties. The length is different and much more 
stable than during the previous V/R active phase in 1955 - 1980. 

Table|2]shows that the cycle length varies and suggests a cor- 
relation between the cycle length and the duration of the triple- 
peak epoch. In fact, the differences in cycle length seem to be 
dominated by the duration of the triple-peak phase. For instance, 
the lengths of Cycles II and III differ by 291 days. But after sub- 
traction of the duration of the triple-peak epochs, the remainders 
of the cycles are of the same length (1033 and 1027 days, respec- 
tively). The lack of an explanation of the variation of the duration 
of triple-peak phases may limit models for the basic V/R activity 



3 For very optically thick disks, the polarization can be parallel to the 
disk plane. 



Stefl et al.: Cyclic disk variations of the Be star f Tau 



15 



as well. Fortunately, at less than ±10 %, the quantitative effect is 
small. 

The well-observed V/R maxima in Ha have been decreasing 
through Cycles I to III. This may be indicative of a decrease also 
of the amplitude of the perturbation causing the V/R variations. 

7.5. V/R phase differences between emission lines 

The relative shifts in the V/R variations of Balmer, Brackett, 
He i, and metal emission lines should probe the disk at de- 
creasing distances from the central star and so place consider- 
able constraints on any explanation of the cyclic V/R cyclic of 
( Tau (and other Be stars). The phase lag between Ha and the 
higher Brackett lines is of the order of A(f> ai 0.25 (see Fig. BJ. 
Unfortunately, any quantitative modeling is somewhat hampered 
by the medium spectral resolution of both the interferometry and 
IR-spectroscopy. 



region was close to the line of sight, whereas the bisector lines 
of the regions of formation of emission lines requiring higher 
proximity to the exciting central star passed the line of sight up 
to one year earlier. 

7. 7. Narrow- and broad line shell absorptions 



From a subset of the present IR spectra, Wisniewski et al. 
( |2007| l derived the opposite V/R phase relation between IR 
lines and Ha. The larger phase and wavelength coverage of the 
present study appears to exclude this possibility. 

7.6. Triple-peak Ha profiles 

The present compilation of spectra maps the evolution of 
the triple-peak profiles in Cycles III and IV in much detail. 
Qualitatively, it follows a very similar pattern in either cycle. 
But the duration of the triple-peak phases is variable and may 
last from 200 to 500 days. They are not phase-locked to the com- 
panion star. However, this does not invalidate the small-number 
statistics produced by Rivin ius et al.| ( |2006[ ), who report that 
triple-peak line profiles only occur in shell stars that are both 
V/R variable and multiple. 

As is apparent from Fig. [3] the triple-peak profile in Ho- is 
accompanied by a somewhat disturbed profile in H/3. But other 
non-hydrogen emission lines in the visual range are largely un- 
affected. Although the lower resolution of the IR spectra leaves 
some space for undetected line profile deformations, IR emis- 
sion lines, too, (including those of H i) do not seem to undergo 
related variations. 

The absence of triple peaks from optically thin emission 
lines and in particular the stability of 1 8446, which probes the 
Ha-forming region due to Ly/3 resonance pumping, suggest that 
the triple-peak feature does not probably correspond to an ac- 
tual density structure. Triple peak profiles have so far only been 
reported in shell stars. This detail, too, would not easily be ex- 
plained by real density enhancements. 

This gives rise to the speculative conjecture about a change 
in the local escape probability due to distortions of the local ve- 
locity field, which of course would have an effect only on op- 
tically thi ck lines, like Ha. Another possibility, guided by the 
finding by Ste fl et al. I d2007b that triple-peak emission lines oc- 
cur in binary Be stars, might be some resonance between the 
orbital motion of the companion stars and the orbital motion of 
gas in the outer disk. 

Ha and H/3 are not the only optically thick emission lines, 
though: It is commonly accepted that the Paschen and Brackett 
lines as well as some strong IR Fe n emission lines form closer 
to the star than the first few lines of the Balmer series. Since only 
Ha and H/3 show some form of triple peak, this may point at the 
outer regions of the disks as the locus of the physical variability. 
Such a region was on the line of sight at the time of the AMBER 
observations. At that time, the bisector line of the Ha emitting 



As described in Sect. |3.3| NLG lines are formed in a broad diver- 
sity of ionic species and over a wide range of ionization and/or 
excitation potentials beginning with O i and Na i and extending 
to Fe in. Conversely, BLG characteristics are exclusively asso- 
ciated with transitions expected from comparatively cool condi- 
tions; examples are Mg i or O i. The simultaneous, let alone long- 
term, presence of both Fe in and Mg i shell absorption is very 
unusual. A search of t he large FLASH, HE ROS, and FEROS 
database of Be spectra ( Rivinius et al.|[2003] > furnished no sec- 
ond example. 

The occurrence of lines with persistent mixed BLG and NLG 
characteristics suggests that the circumstellar disk does not have 
a simple and smoothly varying structure. Rather, NLG and BLG 
components seem to form in two spatially separated regions with 
distinct physical conditions. Simple common-sense considera- 
tions lead to conflicting conclusions about the location of these 
regions relative to the central star. 

On the one hand, their higher excitation would place the 
NLG closer to the star and the lower-excitation BLG in the outer 
disk. On the other hand, the line profiles point to the opposite, 
provided that the disk dynamics is crudely Keplerian: The nar- 
rowness of the NLG, which in some UV lines hardly exceeds 
the thermal width, as well as their lower RV amplitudes suggest, 
then, that the region of formation of NLG lines is relatively far 
from the central star. The broader BLG lines would form in the 
inner disk, where the velocity ranges are larger. 

Apart from their width, NLG and BLG lines also differ in the 
variation of their equivalent width with V/R phase. While BLG 
lines are strongest in the second half of the cycle, NLG lines are 
strongest around phase 0.25. 

The radial-velocity variations of both groups have compara- 
ble amplitude through the V/R-cycle (see Fig. [3]>. However, there 
is a phase difference of about A</> as 0.25 with the NLG trailing 
the BLG. 

8. Conclusions 

( Tau is a Be-shell star. Its spectral type is B2 or earlier one. Both 
interferometry and polarimetry demonstrate that its circumstel- 
lar disk is flat and seen edge-on. The assumption that all Be- 
shell stars are observed edge- and equator-on has been vital for 
the analysis by Riviniu s~et al.| ( |2006[ ) of the fractional critical 
rotation rates of Be stars. 

The plane of the disk has remained stable for decades. 
Warping or tilting, as diagnosed in other Be stars by |Hummel| 
(1998 ), may be present but was not detected. The disk seems 
to have been persistent for about a century. But the strength of 
emission and shell absorption lines has varied on time scales of 
years to decades. 

The photocenter of the Bry line emission lies in the plane of 
the disk but is offset from the continuum source. The same could 
not be diagnosed from the weaker He I 2.06 yum emission line. 
This may be a data-quality issue or result from an alignment, at 
the time of the observations, of the photocenter of the He I line 
with the continuum source. 

The Bry result may also resolve a general ambiguity in the 
interpretation of V//?-asymmetric emission lines. From the spec- 
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tra alone one cannot decide whether the V/R asymmetry is only 
in velocity or also in configuration space. The VLTI observations 
show that the latter is true in £ Tau. 

Like some other, but by far not all, Be stars, £ Tau is a 
binary. The companion remains undetected at optical and IR 
wavelengths. An intensive search in a dense series of multi- 
wavelength observations of a very broad range of electromag- 
netic observables did not furnish any effect of the companion on 
the long-term variability or orientation of the disk. The variabil- 
ity of the circumstellar disk can, therefore, be considered intrin- 
sic to the disk itself and the central B-type star. 

The most prominent spectroscopic signature of the disk ac- 
tivity is the cyclic V/R variability. Its amplitude can for decades 
drop below the level of easy detectability. Pronounced V/R varia- 
tions were resumed about 1992. Since 1996, three complete V/R 
cycles were observed. Their similarity in length and amplitude 
has enabled a detailed phenomenological description of the vari- 
ability. 

The formal mean cycle length is 1405 - 1430 days, with 
cycle-to-cycle variations of less than ±10 %. The basic \og(V/R) 
curve for Ha is smooth and symmetric. However, about the V/R 
minimum a perturbation develops, which has been qualitatively 
the same in all three of the most recent cycles. The Ha emission 
profile develops a comb-like structure at its top, which consists 
of three small peaks separated by equally small depressions. The 
duration of this phase varies from cycle to cycle, and the changes 
may be the main reason of the variation in length of the main 
cycle. In the available observations, triple-peak profiles are re- 
stricted to Ha, with traces showing up in H/3. Therefore, the un- 
derlying perturbation may be confined to the outer disk. 

All emission lines partake in the V/R variability. But there 
are line-specific shifts in phase of up to 25% of the cycle length. 
Ha V/R lags behind Bry, which trails the higher Brackett lines. 
Metallic lines typically behave like the higher Brackett lines. 
This supports the notion that the phase lag increases with de- 
creasing density and increasing distance from the central star. 

The shell absorption lines crudely fall into two different 
groups mainly distinguished by the line widths. Low-ionization 
species and low-excitation lines are broader than the ones involv- 
ing higher energies. Some lines also share the behavior of both 
groups so that each group may predominantly arise from dif- 
ferent locations. If the disk motions are basically Keplerian and 
density and temperature mainly drop with distance from the cen- 
tral star, the grouping provides seemingly contradictory diagnos- 
tics of these locations relative to the central star. Both groups un- 
dergo similar radial-velocity variations but the broad lines lead 
in phase by about one-quarter of a V/R cycle. 

The only variability apparently unrelated to both the V/R cy- 
cle and the binarity is the occasional presence of infrared Ci 
emission. 

Discrete major mass loss events were not observed. But ob- 
servations of other Be stars and the relative constancy of the to- 
tal emission strength suggest that minor ones have likely taken 
place. They may provide a partial explanation of not fully repet- 
itive variations. 

The complexity of the observations and especially the vari- 
ous cyclically repeating phase differences show that the under- 
lying perturbations are not both point-symmetric and radially 
monotonic. Consequently, the perturbations are hardly reconcil- 
able with roughly spherically symmetric models as, e.g., the one 
of Doazan & Thomas] ( |l987| l. Major outbursts might temporarily 
lead to such a situation but the cyclic variability would probably 
require them, too, to be cyclic. A perturbation that propagates in 
a spiral-like pattern has a physical foundation in global m=l disk 



oscillations, which will be developed into a fully self-consistent 
model in Paper II (Carciofi et al., this volume). 
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